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Abstract 

Operator-Schmidt decompositions of the quantum Fourier transform on 
are computed for all Ni,N2 > 2. The decomposition is shown 
to be completely degenerate when A''i is a factor of and when A''i > A^2- 
The first known special case, A^i = = 2", was computed by Nielsen in 
his study of the communication cost of computing the quantum Fourier 
transform of a collection of qubits equally distributed between two par- 
ties. [M. A. Nielsen, PhD Thesis, University of New Mexico (1998), 
Chapter 6, arXiv:quant-ph/0011036 ] More generally, the special case 
Ni = 2"! < 2"^ = was computed by Nielsen et. al. in their study of 
strength measures of quantum operations. [M. A. Nielsen et. al, (accepted 
for publication in Phys. Rev. A); arXiv:quant-ph/0208077 ] Given the 
Schmidt decompositions presented here, it follows that in all cases the bi- 
partite communication cost of exact computation of the quantum Fourier 
transform is maximal. 
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1. Introduction 



Operator-Schmidt decompositions are useful for quantifying the non-local na- 
ture of operators on finite-dimensional bipartite Hilbert spaces. The first special 
cases of Schmidt decompositions of the quantum Fourier transform were com- 
puted by Nielsen to illustrate his study of coherent quantum communication 
complexity. He considered the following problem: 

Suppose Alice is in possession of m qubits, Bob is in possession of 
n qubits, and they wish to perform some general unitary operation 
U which acts on their m + n qubits. How many qubits must be 
communicated between Alice and Bob for them to achieve this goal? 

Nielsen proved that the number Qq (U) of such qubits was bounded by 

1/2 X ifHar (U) < Qo (U) < 2 miu (n, m) , (1) 

where the Hartley strength ii'Har satisfies 

XHar ([/) = l0g2 (Sch(t/)), 

where Sch (U) , defined in Definition^ below, is the number of nonzero Schmidt 
coefficients of U. The upper bound of ^ is trivial, for Alice could simply 
send her qubits to Bob and let him send them back after performing U, or 
vice-versa. To illustrate his theorem, Nielsen considered the quantum Fourier 
transform J-2"x2'^ on n+n qubits. He showed that i^Har (■^2'»x2") = 2n, yielding 
n < Qq {^2^x2") < 2n. Subsequent work by Nielsen |2] improved the general 
lower bound of (pQ) by a factor of two,^ in particular implying that 

Qo (^2"X2") = 2?1. 

In a later paper 13] , Nielsen and collaborators further employ operator Schmidt 
decompositions in the quantitative study of strength measures of the nonlocal 
action of unitary operators.^ Besides revisiting the Hartley strength, among the 
several strength measures considered is the Schmidt strength, 

where U is a unitary operator on TL® IC, {Xk\ are its Schmidt coefficients, and 
H is the Shannon entropy. They give a Schmidt decomposition of J^2"»x2" on 
ni + n qubits for the case m < n and conjecture that Sch (jF2mx2'») = 2^" for 
m > n. 

^ See also footnote |^ for a brief outline of an alternative proof. We remark that Nielsen 
considers qubits for convenience only. In particular, let V be a unitary on C^i ®C^^ , where A^i 
and N2 are the respective dimensions of Alice and Bob's quantum states, with no requirement 
that A^i and N2 be a powers of two. Then any quantum computation of V employing some 
combination of qudit communication and ancillae, possibly of varying dimension, satisfies 
the following bound: X]di2^d'°§2 i<^) ^ ^har where A^^ is the number of qudits of 
dimension d communicated between Alice and Bob. It is assumed that at the end of the 
computation that Alice and Bob retain posession of their (now altered) data qudits, although 
the bound holds whether or not a given net transfer of the (restored) ancillae is allowed. 

^They also consider more general quantum operations than unitaries. 
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1.1. Results 



Schmidt decompositions of the quantum Fourier transform T]^^xN2 ■ 

C^'^ (g) are given for all Ni,N2 > 1, with no requirement that 
either A^i or N2 be a power of two. As a special case, the conjecture of Nielsen 
and collaborators is affirmed. In all cases, the results of Nielsen imply that 
the bipartite communication cost of exact computation of the quantum Fourier 
transform is maximal. Once stated, the decomposition is easily verified; a short 
derivation is given in the Appendix. 

1.2. Definitions and Notation 

Definition 1 Let N, Ni, N2 be integers greater than one satisfying N = NiN2. 
The quantum Fourier transformation^ Tn '■ — > is the unitary 
operator satisfying 

= ^ exp y-^ts^ \t)j, , s e {0, AT - 1} , 

where each |s)^ is a standard basis element. The quantum Fourier transforma- 
tion J^NixN2 on C^^ ^C^^ is obtained by identifying with C^^ (gC^^ under 
the mixed- decimal representation, which asserts the equalities 

when 

s = kN2+e,ke {0, Ni-i},ee {0, N2-1}. 

Remark 2 In the case that Ni ^ N2, the reader is warned that the operator 

J^Ni X N2 equivalent in what follows to Tjss^ x Nx ■ Specifically, !Fn does not 

commute with the unitary operator Rni,N2 ■ given by 

Rn^,N2 \kN2 +£) = \£Ni + k), e {0, Ai - 1} , ^ e {0, 7V2 - 1} , 

which interchanges the digits of the mixed-decimal representation. 

Notation 3 Let H be a finite- dimensional Hilbert space. Then B (Ti.) is the 
Hilbert space of linear transformations on H with the Hilbert- Schmidt inner 
product {A,B)jj(^-^ = TrA+B.* 

Definition 4 Let Ti. and K be finite- dimensional Hilbert spaces, and let F be a 
nonzero linear transformation onH®IC. An (operator) Schmidt decompo- 
sition of F is a decomposition of the form 

Sch(F) 

F= J2 ^kAk^ Bk, Xk > 0, (2) 

k^l 

*This is unitarily equivalent to the discrete Fourier transform. 

''if A is a linear operator on H, then is defined by (x, Ay)j^ = (^A^x,y).j_^ for aXlx,y E Ti,. 
Here (•, •).^ is the inner product on Ji, and we will always take inner products to be linear 
in the second argument. 
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where {^4^}^.^-^ sch(_F) '^"'^ {^fc}fe=i Sch(F) orthonormal sets^ of operators 
on Ti. and )C respectively, under the Hilbert- Schmidt inner product. The quantity 
Sch (F) is called the Schmidt number, and the Xk are called the Schmidt 
coefficients. Such a decomposition is said to be completely degenerate if 
Sch ~ (min (dimTi, dim/C))^ and all the Xk are equal.^ 

We remark that the operator-Schmidt decomposition is just a special case 
of the well-known Schmidt-decomposition 

Sch{ip) 

ip ^ ^ XkCk'S) fk, Afc > 
fe=i 

of a vector ip G TCq (8) JCq, where the {e^} and {/&} are orthonormal7 In par- 
ticular, one sets Ho = B {H) , JCq = B (JC), and = F £ B {H ^ JC). The de- 
composition Q is then obtained by identifying B {Ti) (8) B (/C) with B [Ti,® JC) 
under the natural isomorphism.^ It follows that F and G va B {H ® JC) ~ 
B [Ti) (g) B (JC) have the same operator- Schmidt coefficients, counting multiplic- 
ity, iff 

A= {V®Y)B 

for some unitary "super-operators" V E B {B (H)) and Y E B {B (/C)).^ 

The well-known procedure for computing Schmidt decompositions is re- 
viewed in Theorem |H1 of the Appendix. We content ourselves here with the 
statement that the Schmidt coefficients of V' S Ho (Xi /Co are the square roots of 
the nonzero eigenvalues of the reduced density matrix 

= Tr IV-) (V'l ■ 

Equivalently, the Schmidt coefficients are the nonzero singular values of the 
operator B^ -.Ho-^IC^ given by 

{BiJHd) ^{i^,f<» 9)no®iCo ' 

^but not necessarily bases 

''More generally, if F : H^K — > H' , then one may consider decompositions of the form 
O . where now the Aj^ : H — > H' and : K ^ K' are orthornormal. A useful such decompo- 
sition exists for the communication operator C : (C"i (gi C"2) (g, C"3 C"i ® (C"2 (g, C"3), 
defined by C (a b) (g c = a (b c). One may check that C = \/"T?i3^fc ® B^., 

where = n'^^"^ TJlLx l«> : C"! ® C"2 ^ C"i and Bfc = n'^''^ YJtl\ l^«> («l ■ ^"^ ^ 
(gjC^ , Replacing the swap operators in section III.B.3 of by communication operators, 
one obtains the aforementioned sharp quantum communication complexity bound of 
''^See 13 for a discussion of the Schmidt decomposition. 

*In particular, there exists a unique unitary S : B {H) B (K.) —> B (Ti.® K.) such that 
(S [A®B)) {f ®g) = (Af) ® (Bg) for all / e W and g & K. Here ® denotes the defining 
tensor product of B (H) (g B (K,), considering B {Ti.) and B (fC) as abstract Hilbert spaces. 

^See exercise 2.80 of 1^. One would like to know much more, i.e. invariants which 
specify when are there local unitaries U,Y £ B (Ti.) and V,Z £ B (K.) such that A = 
{U Cgi S (y Z). Such invariants are known only in the two-qubit case, |S| where one 
has the corresponding canonical decomposition of Khaneja, Brockett, and Glaser |S] (see also 
Kraus and Cirac T for simple "magic basis" proof.) 
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where /Cg is the dual space of continuous Uncar functionals on ICq-^'^ 

2. Schmidt Decomposition of 

Notation 5 Let Zn, = {0, iVi - 1}, Zjv^ = {0, iVa - 1}, = Zat, x 
Zpf^, Nil? = {(Nix, Niy) \ x,y ^ Z}, \x~\ = min{n G Z | n > a;}, and [xj = 
— [—a;] . Denote the cardinality of a set C by \C\. Its characteristic function 
Xc satisfies 

, , f 1 if X e C 
^^(^) = \ ^f x^C ■ 

Adopt the convention n modm G Zm- 

Theorem 6 Define an equivalence relation ^ on Z^^ by 

m £- in e NiZ^, 

where the subtraction is not modular, and define M = Z^^/ ~ to be the set of 
equivalence classes}^ Then a Schmidt decomposition of TnixN2 given by 



^wixjv. = E V ^ 1^1 » Be, (3) 
ceM * ^ 

where the matrices of Ac ■ and Be ■ — > C"'^^ are defined by 

(^c)feife2 = TT^exp ^ {N2kik2 + kiC2 + k2Ci) ki,k2eZN, 

with each (ci,C2) G C arbitrarily chosen. (Ac doesn't depend on this choice.) 

^"See 13 for a proof that the Schmidt decomposition is a consequence of the singular value 
decomposition. In fact they are mathematically equivalent. 

^^The reader may check that for A^i = 2 and N2 = 3 that Ai consists of 
{(0,0), (0,2), (2,0), (2, 2)}, {(1,0), (1,2)}, {(0, 1) , (2, 1)}, and {(1,1)}. 
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Proof. It is trivial to check that {Ac} and {Be} arc orthonormal sets. 
Furthermore, for ki,k2 G ZjVi and £±,£2 € Z^r^, 

{kiJil ( E JW\C\Ac®Bc] \k2,£2) 

= E y^iq (fell Ac |fc2)(^l| i?C K2) 



y — 



E 

1 



exp 



exp 



exp 



27ri 



2TTi 

— (iV|fci/c2 + iV2fcl£2 + N2k2Ci + ii£2) 



2-Ki 



{Nikik2 + N2ki£2 + N2k2£i + 4^2 



xcm,£2)) 

xc{{h,£2)) 



{N2ki + £1) {N2k2 + I2) 



= {ki£i\TNiy.N2 1^2,^2) , 

as desired. ■ 

The reader may find it instructive to compute the linear spans of the matrices 
Be corresponding to each of the Schmidt coefficients. 

Corollary 7 The Schmidt decompositions of J-'nixN2 /^^^^ three categories: 



1. If Ni is a factor of N2, then there is only one Schmidt coefficient, \J N2 /Ni , 
with multiplicity . 



2. If Ni > N2 , there is only one Schmidt coefficient, yj N\ /N2 , with multi- 
plicity N2 ■ 

3. Otherwise, J-N1XN2 has three distinct nonzero Schmidt coefficients: 



N2 


Ni 




Ni 


N2 




~ N2 




N2 


Ni 


Ni 




_Ni_ 


N2 


N2 


2 

Ni 




_Ni_ 


N2 





, multiplicity mod A^i)^ 



, multiplicity ((— A^2) mod iVi 



In all cases, the Schmidt number of J^NixN2 min (iV^, iV|) . In particular, the 
Schmidt decomposition is completely degenerate in Cases 1 and 2. 

We remark that the previously known cases fall under Case I. Case 2 verifies 
the Schmidt numbers conjectured in j^. Since the Schmidt decomposition in 
Case I (or Case 2) is completely degenerate, Theorem|Hl (below), may be used to 
find a Schmidt decomposition of the form of equation (0) for any orthonormal 
basis {Ak} (or {Bk}, in case 2)}"^ 

^^Note that cases 1 and 2 overlap for A^i = N2- 
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3. Appendix: A Derivation 

It will soon be apparent that the crucial fact which allows easy calculation of a 
Schmidt decomposition of T is the following: No two of the Be have a nonzero 
matrix entry in the same place. 

The well-known computational recipe needed here is summarized in 

Theorem 8 Let t/j £ H ^ K, be nonzero. If 

PK^Trli,) {i;\=J2f,,\f,) {f,\, 
rt — 

is a spectral decomposition of the reduced density matrix, then a Schmidt de- 
composition of is given by 

{^lw>o} 

where each ei is defined by the requirement that 

= (5) 

for all V Cz Ti. Furthermore, all Schmidt decompositions of ^ may be exhibited 
in this manner. 

Derivation of Theorem |6l We follow the prescription of Theorem |H1 
and employ the natural isomorphism B (C^i) (g) B (C^^) ~ B (C^i (g) C^^)^ as 
explained in section [l.2.I The reduced density superoperator p £ B (^B (C^^)) is 
defined by the equation 

E 

for arbitrary A,B £ S(C^"), where E runs over a basis of B (C^i). For 
i E and i E define the standard basis elements 

Ej = \n) O2I e B (C^O ' 1^1) (^2| e B (C^^) . 

We compute p by studying its matrix coordinates 

Pern = (^f'/'-^™)s(c"2) ' 
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Similarly, let 
Then 

Pirn ~ -^Jr^Jm 



J 62 



^ ^1 exp + £i) {N^n + ^2)) \ 

TV ^ I xexp -^(iV2ji+mi)(7V2j2 + m2)) ^ 



exp I — [liit2 - mim2) 



X XI *^^P ( - ^2) ii j X exp ( ^ (^1 - "^i) J2 

ii=o ^ ^ ^ j2=o ^ ^ 



Evaluating the appropriate inverse-Fourier transforms, 



(^1^2 - rmma) j x xn,z^ mj . (6) 



The spectral decomposition of p into a linear combination of projections may 
be simply read off from the asymptotic n 00 behavior of © to the power of 
n E Z+.^^ One need not do this, however, for using the identity 



XNi 



/iZ^ (i- mj = Y xc XC {rn) 
equation © may be rewritten as 



^= E ^\C\x\Bc){Bc\, 
ceM 2 



where the Be are orthonormal, as noted before. The Ac are easily computed 
using ■ 
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